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Strong transmission and reflection of edge modes in bounded photonic graphene 

Mark J. Ablowitz^, Yi-Ping Ma^Q 
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The propagation of linear and nonlinear edge modes in bounded photonic honeycomb lattices 
formed by an array of rapidly varying helical waveguides is studied. These edge modes are found 
to exhibit strong transmission (rehection) around sharp corners when the dispersion relation is 
topologically nontrivial (trivial), and can also remain stationary. An asymptotic theory is developed 
that establishes the presence (absence) of edge states on all four sides, including in particular 
armchair edge states, in the topologically nontrivial (trivial) case. In the presence of topological 
protection, nonlinear edge solitons can persist over very long distances. 

PACS numbers: 42.70.Qs, 42.65.Tg, 05.45.Yv 


There is significant interest in using optical materi¬ 
als to realize phenomena that are difficult to observe in 
traditional materials. In particular, light propagation in 
photonic lattices with a honeycomb (HC) background, 
referred to as photonic graphene due to its similarities to 
material graphene m , has been extensively studied both 
experimentally |2|, [3| and theoretically [j, l5[. In infinite, 
or bulk, lattices the wave dynamics exhibit the interest¬ 
ing feature of conical diffraction [2|. This phenomenon 
has been explained by the presence of Dirac points, or 
conical intersections between dispersion bands 

Recently it was shown that introducing edges and suit¬ 
able waveguides in the direction of propagation, unidirec¬ 
tional edge wave propagation occurs at optical frequen¬ 
cies 0 . The system is described analytically by the nor¬ 
malized lattice nonlinear Schrodinger (NLS) equation 

idzfj = -(V + -h V{r)'ip - (Tq \ fj, (1) 

where the scalar field is the complex envelope of the 
electric field, z is the direction of propagation and takes 
on the role of time, r = (x, y) is the transverse plane, V = 
{dx^dy), the potential V{t) is taken to be of HC type, 
and the coefficient <Jo is the strength of the nonlinear 
change in the index of refraction. The vector field A{z) 
is determined by helical variation of the HC lattice in the 
direction of propagation, and plays the role of a pseudo- 
magnetic field. The particular choice used in [7| is, in 
terms of dimensionless coordinates, 

A(z) =/^(sinf^z, — cosf^z), (2) 

where k and Ul are constant. The waveguides defined by 
A{z) are written into the optical lattice using a femtosec¬ 
ond laser writing technique [H. The associated linear 
edge wave propagation was investigated experimentally 
and computationally in the tight binding limit in Q . Re¬ 
markably, these edge waves were found to be nearly im¬ 
mune to backscattering. This phenomenon was found 
to be related to symmetry breaking perturbations which 
separate the Dirac points and leave a nontrivial integer 
“topological” charge on the separated dispersion bands. 
In this setting, photonic graphene exhibits the hallmarks 
of Floquet topological insulators @ . 


Motivated by this work, the existence of traveling 
edge modes was investigated analytically in [l^ for gen¬ 
eral periodic pseudo-fields A{z) and with nonlinearity, 
i.e. (Jo 7^ 0 in the tight binding limit. Assuming that 
the pseudo-field varies relatively rapidly, which is consis¬ 
tent with the experiments in [7|, an asymptotic theory 
based on Floquet theory was developed which leads to a 
detailed description of the wave dynamics. Importantly, 
in the presence of weak nonlinearity the classical one¬ 
dimensional NLS equation was found to govern the enve¬ 
lope of traveling edge modes. When the NLS equation is 
focusing a family of nonlinear edge solitons exist due to 
the balance between dispersion and nonlinearity. When 
the dispersion relation is topologically nontrivial, these 
nonlinear edge solitons are immune to backscattering as 
with linear traveling edge modes. The existence of these 
nonlinear edge solitons vastly expands the landscape and 
understanding of localized states along edges. 

Much of the existing literature on edge modes in HC 
lattices, including the above referenced work, has focused 
on extended edges of either zig-zag or bearded type. Here 
we consider a rectangular-type domain, which has zig-zag 
edges on the left and right, and armchair edges on the 
top and bottom, as illustrated in Fig.[Tl In the absence of 
helical waveguides, armchair edge states do not exist in 
isotropic HC lattices ini- However, with rapidly varying 
helical waveguides we find traveling edge modes can exist 
along the armchair edges. There are two types of arm¬ 
chair edge states corresponding respectively to isotropic 
and anisotropic HC lattices. In the presence (absence) of 
these armchair edge states, linear edge modes and nonlin¬ 
ear edge solitons are found to exhibit strong transmission 
(reflection) around the sharp corners under certain con¬ 
ditions; otherwise they may decay due to dispersion or 
scattering. The ability to easily switch between trans¬ 
mission and reflection as well as remaining stationary, by 
varying the underlying parameters provides a new means 
for the control of light conferred by appropriately merg¬ 
ing topological, linear and nonlinear effects. 

To employ the tight binding approximation, we sub¬ 
stitute into Eq. ([T]), and assume a Bloch 
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FIG. 1: A 5 X 5 bounded HC lattice with zig-zag edges 
on the left and right, and armchair edges on the top 
and bottom. The A and B sites are indexed following 
Ref. [HI. The primitive lattice vectors vi and V 2 and 
the intersite distances dj, j = 0,1, 2 are also labeled. 


wave envelope of the form Q 

(j) ^ ^ ^ (<^v(^)01,v T ^v('^)02,v) (3) 

V 

where = (/)i(r — v), 02,v = 02(i* — v) are the lin¬ 
early independent orbitals associated with the two sites 
A and B where the HC potential R(r) has minima in 
each fundamental cell. The lattice sites are located at 
V = /ivi -i- / 2 V 2 , where I 1 J 2 ^ and the lattice vec¬ 
tors vi and V 2 are given by vi = (\/3/2, 1/2) , V 2 = 
(v^/2, — 1 / 2 ) . Carrying out the requisite calculations 
(see [3] for more details), we arrive at the following two 
dimensional discrete system 

idzCimn T (^—(' 2 ^)^)mn T ^mn — O 5 ("1) 

id zb 

mn + ('^+{^)«) mn ^mn — O 5 (5) 

where a depends on the underlying orbitals, is propor¬ 
tional to (Jo, and the linear operators jC± take the form 

{C-h)mn ^m+l,n-l )5 

(/^_|_ < 2)77271 ^ ^mn T pi^ ^ ^772+1,77+1 T 6 ^ ^772 —1^72+1 ) 5 

where p is a lattice deformation parameter, and Oj = 
dj • A, j = 0,1, 2. The vector do is the direction vector 
from the site fo the site and is approximately 

{l/\/3,0) in the tight binding limit US- The vectors 
di ,2 = do — vi ^2 are the direction vectors from any given 
B site to the other two neighboring A sites. For later 
reference we express the full 2D discrete system (HH) in 
the compact form 

'^dzCjYiji {/^cjrnn T <^ 1 ^ 77777 ,! C 77272 = 0, ( 6 ) 


where 


C = 


0 C-{z) 

C^{z) 0 


The HC lattice is formed by those sites with m-\-n even 
as shown in Fig. [TJ 


Using the pseudo-field given by Eq. ([2]) with period 
T = 27r/U = 1.5 (consistent with the experiments in Q), 
our 2D discrete system ([6]) has two free parameters p 
and cr; p = 1 and cr = 0 for a perfect HC lattice without 
nonlinearity in Q. The initial condition is taken to be 
localized on the left zig-zag edge of the form 

Cmn{z = 0) = C{m, z = (7) 

where ujq is the carrier wavenumber in m (along the 
edge), and is the ID edge state that decays in n (per¬ 
pendicular to the edge). The envelope function C, whose 
evolution is determined below, is taken to be 

C(m, 2 ; = 0) = Asech(z/(m — mo)), (8) 

where A is the amplitude, v is the spectral width , 0 < 
z/ ^ 1, and mo = 0 unless otherwise stated. The 2D 
evolution can be visualized using space-time plots of the 
amplitude \b{z)\ along the four edges, as shown in Fig.O 
The effective mass of the mode for any z is the square 
of the L^-norm in an interval of width 4/z/ around the 
center of mass evaluated at the edge. 

First we consider the linear case (a = 0), and choose 
the initial condition to have c^o = 7r/2, u = 0.2, and 
A = 1 consistent with Q- In this case we find that two 
contrasting types of linear wave dynamics can exist in 
the bounded HC lattice depending on the choice of p. As 
shown in Fig. EJ^a), for p = 1 > 1/2, the edge mode is al¬ 
most perfectly transmitted around the sharp corners and 
thus loops around the HC lattice. The effective mass re¬ 
maining after four (eight) loops, or 2 ; = 1550 {z = 3100), 
is 96% (94%). As shown in Fig. Efb), for p = 0.4 <1/2, 
the edge mode is very strongly reflected at the sharp 
corners and thus bounces back and forth along the left 
edge. The effective mass remaining after two (four) re¬ 
flections, or 2 ; = 1250 {z = 2500), is 96% (89%). See 
[HI and Q foi* animations of Fig. EJ^a) and Fig. El^b). 
For either choice of p, choosing cjo % 7r/2 causes the edge 
mode to disperse more and thus destroys coherence. Nu¬ 
merical results not shown here indicate that the overall 
wave dynamics are insensitive to the detailed shapes of 
the corners, though the edge mode lingers longer at the 
corners containing degree-1 vertices (e.g. the upper and 
lower left corners in Fig. [T]). 

To describe a zig-zag (armchair) edge mode we take a 
discrete Fourier transform in m (n) on Eq. (|6|), i.e. letting 
Cmn = {Cmn = Eollowing 0, llS ^ WC 

assume that the pseudo-field A or equivalently the linear 
operator jC varies rapidly, namely jC = >C(C), where ( = 
e~^z^ |e| << 1, and thus has period T = e~^T. Then 
we can employ a multiple-scales ansatz c = c(z,(), and 
expand c in powers of e as c = + * * + where 

c denotes the vector Cmn, Cn oi* Cm depending on the 
context. At 0(e“^), = 0 which leads to = 

€^^^( 2 ;). At 0(1), to remove secularities, we get 

idzC^^^ + = 0 , 


(9) 
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(a) (b) 


FIG. 2: Transmission and reflection of linear edge 
modes (a = 0) in a 65 x 65 bounded HC lattice, shown 
as space-time plots of \b{z)\ at the edges. Here ‘zz’ 
represents zig-zag and ‘ac’ represents armchair; the 
horizontal axis begins with the left zig-zag edge with 
increasing m. The pseudo-field is given by Eq. @ with 
K = lA and period T = 1.5, and the deformation 
parameter p is (a) p = 1; (b) p = 0.4. The initial 
condition 03 has ujq = tiI2^ v — 0.2, and A = 1. 

See 0 and [ij for animations of panels (a) and (b). 

where £ is given by the Magnus-type expansion [lH 

£ = £^^^ + e£^^^ -h + O(e^), (10) 

with m = f-^ £ c{0dc, 

= A r [\c{c)C{C) - c{C)C{0)dCdc 

2i Jo Jo 

and £^‘^^ can be written in terms of a triple integral not 
explicitly shown. The discrete system m conserves both 
the mass N and energy H given by 

7V= (c(0),c(0)), (11) 

H = - (c(°),(t|c(°)|2c(°))/2, (12) 

where the inner product is defined as (f, g) = fj9j- 
Assuming that the spectral envelope is centered at ujq 
with width z/, 0 < z/ ^ 1, the edge mode may be ex¬ 
pressed in terms of a scalar envelope function C(z;cj) as 

= C(z; cj)c^ (13) 

where is the ID edge state given by 

£(u;)c^ = -a(oj)c^, (14) 

normalized such that ||c ^||2 = = 1- The func¬ 

tion a{(jj) is the dispersion relation. The solvability con¬ 
dition on Eq. ([9|) leads to 

id^C - a(uj)C + a{uj)\CfC = 0, (15) 

where a{uj) = (Tani{uj) with ani{uj) = Hc^Hl = T* 

When the nonlinearity a ^ 0 and second order dis¬ 
persion a"{ujo) 7 ^ 0, maximal balance between dispersion 
and nonlinearity leads to 

^0^0 + olqLU + C dolCpC = 0, (16) 


where uj = uj — c^o, (co’o), J = 0,1,2, and 

do = d(co’o)- The group velocity of the envelope is a'^. We 
remark that the special stationary case can occur when 
0^0 = 0; in general if the pseudo-field A is independent of 
z, then the time-independent linear operator £ can admit 
the so-called “zero-energy” edge states with a{uj) = 0 
identically 0. Replacing uj by —id^ where y repre¬ 
sents the m (n) direction for the zig-zag (armchair) edge, 
Eq. (pTj) becomes the ID classical second-order focusing 
(defocusing) NLS equation when a'^cr > ^ (<^o^ < 0). In 
the focusing case, Eq. m admits a two-parameter fam¬ 
ily of solitons [l6| ; the two parameters may be chosen as 
ujQ and z/. A 2D edge soliton is then obtained from the 
initial conditions ([7]-[8|) with 

A = vyJa'Q/do- (17) 

Thus robust edge solitons exist for a finite interval of c^o; 
at each cjo, the choice of v fixes the amplitude A. 

The zig-zag (armchair) dispersion relation can be nu¬ 
merically computed on a ID domain with suitable bound¬ 
ary conditions on both ends and compared with analytic 
predictions based on Eq. m- Eigure [3] shows the dis¬ 
persion relations computed using the same parameters as 
in Eig. [2] and plotted on the periodic interval uj G [0,7r). 
The Eloquet eigenvalues corresponding to the eigenfunc¬ 
tions with at least 90% of the mass concentrated on the 
left (right) half domain are highlighted in red (green). 
The Z 2 topological index / is defined as the number of 
intersections modulo 2 between either the red or green 
curve and the horizontal axis 0 . Interestingly, in some 
cases there exist edge states near the outer boundaries 
of the bulk spectrum at large |(a|, but since they are not 
related to the wave dynamics at small \a\, they will not 
be discussed further here. 

As established in 0 for the zig-zag edge, to leading 
order the ID edge states satisfy 

£(0)c^ = 0. (18) 

Eor those edge states localized on the left (right), the 
mass is concentrated in the b (a) sites. The decay expo¬ 
nents in n are 0(1), with amplitudes 0(1) on the edge. 
The dispersion relation a{uj) is 0(e) and is given by 

a{uj) = (19) 

thus the group velocity of the envelope is 0(e). The 
zig-zag dispersion relation computed using the same pa¬ 
rameters as in Eig. [2](a) (Eig. [2l^b)) is shown in Eig. [31(a) 
(Eig. |3l(c)). This dispersion relation is topologically non¬ 
trivial (/ = 1) for p = 1 > 1/2 and topologically trivial 
(/ = 0) for p = 0.4 < 1/2. In the former case, the 
left zig-zag dispersion relation terminates at two points 
denoted respectively by (cc;^_,a_) near the upper bulk 
and {uJm-\-^ fhe lower bulk, while the right zig¬ 

zag dispersion relation has the opposite sign. These four 
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FIG. 3: Dispersion relations on the (a,c) zig-zag and 
(b,d) armchair edges. Values of A and p in panels (a,b) 
and (c,d) agree respectively with Fig.[2fa) and Fig.[2fb). 


points, known as the massive Dirac points, result from 
splitting the massless Dirac points in 

For the armchair edge, describes isotropic 

graphene when A is invariant under 27r/3 rotations, 
and anisotropic graphene otherwise. Thus the circular 
pseudo-field given by Eq. ([2]) leads to isotropic graphene. 
Anisotropic graphene may be exemplified using the ellip¬ 
tic pseudo-fields invariant under tt rotations introduced 
in [l0|. In this case the edge modes behave similarly to 
those in the zig-zag case. 

In the isotropic armchair case, does not admit 
edge states [uj. However, as shown in Fig. [SJb) and 
Fig. md) computed using the same parameters as in 
Fig. EJa) and Fig. EJb), armchair edge states exist near 
the inner boundaries of the bulk spectrum. The topo¬ 
logical index of the armchair dispersion relation is found 
to be the same as the zig-zag case, namely / = 1 for 
p > 1/2 and / = 0forp<l/2. The existence of these 
edge states in L can be explained by expanding C to 
O(e^) in Eq. ([T^ . Eor these armchair edge states, the 
mass is equally distributed in the a and b sites. The de¬ 
cay exponents in m are 0{e)^ with amplitudes 0(e^/^) 
on the edge. The dispersion relation a{uj) is 0(1), and 
so the group velocity of the envelope is also 0(1). The 
detailed calculations will be presented in a separate pa¬ 
per. Here we note that in the 1 = 1 case as shown in 
Eig. mb), the armchair dispersion relation does not ter¬ 
minate at the massive Dirac points, which are located at 
(0, ±a_) and (0, ±(a+) and thus are degenerate in 

The linear and nonlinear dynamics of edge waves can 
be understood in terms of the dispersion relations and the 
conservation of mass N and energy H defined in Eqs. (pTF 


[E]). Consider a J-soliton configuration where each soli- 
ton is given by Eq. © with envelope function Cj , carrier 
wavenumber ooj and carrier frequency j = 1, • • • , J. 
Then to leading order, the mass and energy evaluate to 

N = J2\\Cj\\l H =-Y,c^j\\Cj\\l (20) 

i=i i=i 

where || • ||2 denotes the L^-norm. 

In the linear case, the amplitude \C\ and frequency a 
of each individual wavenumber oo is conserved. Let la 
denote the interval in a spanned by the left zig-zag dis¬ 
persion relation. When p > 1/2, both the zig-zag and the 
armchair dispersion relations are topologically nontrivial 
as shown in Eigs. m^“t)). In this case, corresponding to 
any a e la there is a single value of oo and thus a single 
unidirectional edge mode on each of the four edges. The 
strong transmission of linear edge modes in Eig. EJ^) is 
then explained by a global counterclockwise mode formed 
by ‘gluing’ together these four edge modes through the 
four corners. When p < 1/2, both the zig-zag and the 
armchair dispersion relations are topologically trivial as 
shown in Eigs. m^-d). In this case, corresponding to 
any a £ Xa there are two values of uj and thus a pair 
of counter-propagating edge modes on either of the two 
zig-zag edges. The strong reflection of linear edge modes 
in Eig. EJb) is then explained by the pair of edge modes 
on the left zig-zag edge transforming into each other via 
scattering at the upper and lower left corners. In short, 
the transmission channel is available and the reflection 
channel is unavailable when the dispersion relation is 
topologically nontrivial, and vice versa. We remark that 
for other parameter choices, the transmission and reflec¬ 
tion channels may be simultaneously available, as well 
as the scattering channel into the bulk. These possibili¬ 
ties results in weaker transmission/reflection, and so are 
outside the scope of this paper. 

Next we consider a nonlinear case: a = 0.02 ^ 0, with 
initial conditions c^o = 57r/8, v = 0.2, and A given by 
Eq. (pT|) . In this case the ID NLS equation ([T6|) is focus¬ 
ing and so admits edge solitons on the zig-zag edges. In 
order that the zig-zag dispersion relation spans a wider 
interval in ujm and thus has greater curvature a", we 
choose smaller p in the topologically nontrivial case fol¬ 
lowing 0- As shown in Eig. [41(a) for p = 0.6 > 1/2, 
the edge soliton is almost perfectly transmitted around 
the sharp corners. The effective mass remaining after 
two (four) loops, OT z = 2000 {z = 4000), is 93% (93%). 
As shown in Eig. IlKb) for p = 0.4 <1/2, the edge soli- 
ton is strongly reflected, but noticeably not as well as 
in the linear case, at the sharp corners. The effective 
mass remaining after two (four) reflections, oi z = 1650 
[z = 3300), is 72% (71%). 

In a simple scattering process where an incoming edge 
soliton is scattered by one or more corners into an out¬ 
going edge soliton, the conservation of the carrier fre- 
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(c) (d) 


FIG. 4: Transmission and reflection of nonlinear edge 
solitons (a = 0.02) in a 65 x 65 bounded HC lattice, 
shown as space-time plots of \b{z)\ at the edges. Values 
of A and p in panels (a,c) and (b,d) agree respectively 
with Fig. Ha) and Fig. [Sfb), except that p = 0.6 in 
panels (a,c). The initial condition CZHSl) has A given by 
Eq. ([T7|) . ly = 0.2 and (a,b) ujq = 57r/8. Panels (c,d) 
have two envelopes with (mo,ci;o)- (c) ( —16,27r/3) & 
(16,77r/12); (d) (-16,37r/4) & (16,57r/8). 


quency a determines the carrier wavenumber uj as in the 
linear case, and the conservation of the envelope mass 
||(7||2, which amounts to the conservation of us¬ 

ing Eqs. (jH]) and (pTj) . determines the spectral width v. 
In the transmission case (p > 1/2), since an incoming 
zig-zag edge soliton is scattered by an even number of 
corners into an outgoing zig-zag edge soliton almost per¬ 
fectly, the conservation laws imply that the zig-zag edge 
soliton also maintains its shape almost perfectly. In the 
reflection case (p < 1/2), the zig-zag edge soliton does 
not maintain its shape as well due to multi-scattering at 
the corners and backscattering along the edge. 

Einally we consider soliton interactions in the nonlinear 
case, with the initial condition chosen to contain two zig¬ 
zag edge solitons with different carrier wavenumbers c^o 
and thus different group velocities. In the transmission 
(reflection) case, as shown in Eig. [4](c) (Eig. [4|(d)), the 
two edge solitons are transmitted (reflected) almost inde¬ 
pendently. The collision between these two edge solitons 
exhibits a phase shift, reminiscent of soliton collisions in 
the classical ID NFS equation. 

In this paper, strong transmission and reflection of 
linear and nonlinear edge modes is demonstrated in 
bounded photonic graphene formed by rapidly varying 
periodic waveguides. The transmission (reflection) phe¬ 
nomenon depends on the presence (absence) of topolog¬ 


ical protection on the zig-zag and armchair edges. Eor 
isotropic graphene the armchair edge states have slowly 
decaying profiles. These unconventional edge states pro¬ 
vide an alternative means to localize and transport light 
on armchair edges without additional modifications [l7| . 
Simple scattering of edge modes by sharp corners is ex¬ 
plained via conservation of mass and energy. In the 
presence of nonlinearity, the existence of nonlinear edge 
solitons greatly broadens the landscape of coherent edge 
modes. The synergy between topology and nonlinearity 
makes bounded photonic graphene an ideal candidate for 
robust routing of electromagnetic energy, which is a topic 
of considerable recent interest [HI . 
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